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Abstract. 

Urbanization has become a hot issue in context of environmental and socio-political 
scenarios which is addressed at every forum internationally. The mega cities are considered 
the main origin of socio-economic development which caused to emerge a number of issues 
like biodiversity, environmental degradation, resource consumption, implementation of law 
and order and provision of basic facilities to the general public. The area under investigation 
consists of Lahore, Gujranwala and Sheikhupura. The study site was bounded by 73-75 E 
longitudes and 31-33 N latitudes. We used Landsat satellite data to map Spatio-temporal 
variations in urban sprawl from 1990 to 2019 with a temporal window of 15 years. The Landsat 
data is free, highly reliable and considered as primary source. The classification results show 
that the total area of study site was site was 29355 km’ including 21933km’ were green index 
4595 km’ was under human settlements and 2827 km? was the waterbody in 1990. The 
classification of Landsat image of the year 2005 describes that area of human settlements was 
increased to 9366 km’, the volume of water body was reduced to 2111km° and the vegetation 
was also degraded to 17878km*. Again, the urban area was computed using satellite imagery 
for the year 2019 which was 16105km/? in 2019. Kappa stat proved the accuracy of supervised 
classification what was around 87%. Remotely sensed datasets proved the reliability of Landsat 
satellite images for estimation of urban sprawl during last three decades. 

Keywords: socio-political development, Urban Sprawl, Landsat, Biodiversity, Classification. 
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Introduction. 


Urbanization has become a hot issue in context of environmental and socio-political 
scenarios which is addressed at every forum internationally. About 15% of the world’s 
population was urbanized in 1990 but now it is expected that half of world’s population will 
be urbanized by 2030 [1]. This urbanization growth has devasted urban balance and effected 
the sustainable development. The situation has become worse which lead to severe 
environmental issues and the climate change [2]. This climate change has affected the cropping 


system [3] glacial melting [4], severe flooding and droughts [5]. 


Many studies have indicated a big difference in urban temperature in comparison to 
its outskirts [2,6,7,8]. Removal of vegetation for construction of buildings, malls and plazas, is 
considered the main cause of rise in urban temperatures [9]. Urban transport adds a variety of 
pollutants in the atmosphere through tail pipe emissions which cause to block the reflected 
solar radiations, and finally result into urban heat island [2]. The removal of vegetal cover from 
urban areas cause the solar radiations to approach earth’s surface directly which result to an 
abrupt increase in urban temperatures in comparison to its surroundings [2]. 

The mega cities are considered the main origin of socio-economic development which 
caused to emerge a number of issues like biodiversity, environmental degradation, resource 
consumption, implementation of law and order and provision of basic facilities to the general 
public [10,11,12,13]. These issues need to get resolved on emergency grounds. 

Sustainability and urban development have been widely focused in various studies 
[11,12,14,15]. Some studies have proved that regional sustainability is the key parameter to 
evaluate the credibility and authenticity of urbanization [12,14,15] however, many others believe 
that urban sustainability is oxymoron [11]. Urban sustainability has two main pillars that include 
ecological and social sustainability. Ecological sustainability is concerned with ubiquitous urban 
sprawl while social sustainability relates to urban population growth. Recent urban growth in 
Asian countries put huge pressure on its ecological system which caused to emerge new 
settlements e.g., highways that eaten up agricultural fertile lands, wetlands, forests and other 
lands. The socio-economic sustainability is concerned with two main themes 1.e., economic 
balance and social equity [16,17]. Social inequity is concerned with all the issues associated with 
society that includes social injustice, welfare of migrants [18], rural-urban settlement gap, health 
issues, ait pollution, traffic injury and occupational hazards [19,20,21]. 

In the very first step, urban instability hit the society by disturbing literary, 
mobilization, income, education, political participation and health. Next, these issues become 
prominent when the urbanization effects the carrying capacity of a particular urban society. 

Rapid urbanization and industrialization were faced by various developed countries 
[22]. Urbanization of ancient times, was associated with employment opportunities therefore, 
people rushed and settled themselves in surroundings of urban societies for their livelihood. 
The socio-economic development results in peace and prosperity of a region. 

The main objective of this research was to investigate the urban sprawl in mega cities 
of Punjab Pakistan in context of socio-political issues using Remote Sensing and Geographic 
Information System (GIS). 
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Material and Methods. 
Study site. 

We selected three adjacent districts in Punjab Pakistan where high rates of urban 
sprawl were reported in recent three decades. The area under investigation consists of Lahore, 
Gujranwala and Sheikhupura. The study site was bounded by 73-75 E longitudes and 31-33 N 
latitudes. The study site is mapped in Figure 1. 
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Figure 1. Study site. 
The methodology used in this research is mapped in Figure 2, which describe the 
complete mechanism to compute urban growth from satellite image acquisition to final display 
of thematic maps. 
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LANDSAT Image Acquisition 
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Figure 2. Methodology. 
We used Landsat satellite data to map Spatio-temporal variations in urban sprawl from 
1990 to 2019 with a temporal window of 15 years. The Landsat data is free, highly reliable and 
considered as primary source. We obtained satellite images of Landsat 5 for the years 1990 and 
2005 and Landsat 8 for the year 2019. The complete description of satellite images used to 
accomplish this research is mentioned in Table 1. 
Table 1. Landsat 5, 8 image acquisition dates. 


Sr No. Satellite Date Resolution Bands 
1 Landsat 5 Nov 11, 1990 30 m’ 7 
2 Landsat 5 March 16, 2005 30 m’ 7 
3 Landsat 8 April 01, 2019 30 m’ 11 
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A Landsat image is comprised of a number of spectral bands, commonly known as 
layers. In the first step, these layers were combined to obtain a satellite image using layer stack 
utility in Erdas Imagine 14. This utility is capable to bind these layers carefully without merging 


their information at cell level. An overview of layer stacking is shown in the Figure 3, 
Layer 1 


Layer 2 
Layer 3 


Layer 4 


Figure 3. Layer staking in Erdas Imagine. 

The swath width of Landsat image is 185km^. It is a very large image in comparison to 
the spatial extent of the study site. Therefore, it was important to figure out the investigation 
site from a large dataset for saving the processing time. We used the mask of study site and 
applied “Extract by Mask” to obtain the investigation site only from the Landsat images. 

All satellite images were checked geometrically to cross validate the spatial locations 
of various landuse features. To do so, we marked some features on satellite images and checked 
their locations through Global Positioning System (GPS) but we could not find considerable 
variations. 

Classification is used to demarcate and discriminate between various landuse features 
existing in the study site. The main objective of classification was to compute the volume in 
term of area of existing features. There are two types of classifications, 1) Supervised 2) 
Unsupervised. In supervised classification, we do recognize landuse features to computer by 
selecting trainee samples. Software intake these sample sites and recognize the existing landuse. 
Supervises classicization is considered technically correct in comparison to unsupervised 
classification. We computed the variations in area of landuse features through supervised 
classification. 

Kappa coefficient is considered a good indicator to determine the accuracy of 
supervised classification [23]. It is widely used therefore, we applied Kappa stat for accuracy 
assessment of supervised classification. 

Result and discussion. 
The supervised classification of Landsat images is mapped in figure 4. Figure 4 


describe that urban boundaries have sprawled beyond the limits during last three decades 
which has eaten up fertile agriculture lands. Figure 4(A) is showing the spatial extent of water 
body, vegetation and the built-up area in the year 1990. According to figure 4 (A), the total 
area of study site was 29355 km’ including 21933km° were green index 4595 km* was under 
human settlements and 2827 km’ was the waterbody in 1990. 

Figure 1(B) is showing that the area of human settlements had been increased in 2005 
in comparison to 1990 from 4595km? to 9366 km’, the volume of water body was reduced to 
2111km/* and the vegetation was also degraded to 17878km‘’. 
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Figure 4. Spatio-temporal variations in landuse of the investigation site from 

1990-2019 having a temporal window of 15 years. 

Figure 1(c) was mapped using Landsat 8 image of year 2019 which show a reasonable 
jump in built up area that reached to 16105km* and hence a decline to vegetation as the area 
under vegetation was observed as 11888km/’. Figure 5 is showing temporal variations in built 
up, vegetation and the water body. 
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The accuracy assessment of supervised classification for the years 1990, 2005 and 
2019, we applied Kappa stat on the classified images and applied 25 sample points including 
5 to waterbody, 10 to build up and the vegetation. The calculations are as follows, 


Table 2: Sample points applied to classified image for the year 2019. 


Class Vegetation Build-up Water Total-User 
Other-vegetation 9 1 10 
Build-up Area 10 10 
Water 1 4 5 
Total-Producer 10 11 4 25 


TS= Total number of Sample 
TSC= Total number of corrected Sample 


_ (TS * TCS) — (col. total x Row total) 
Kappa coefficients(K ) = — (TS)? — ¥(col. total * Row total) _ x 100 


ae — [25+ 23) - {10 +10) + (10*1D + 6] aoo 
~~ (25)2—{(10* 10) +(10* 11) +(5*4)} 


_ (S75) — {(100) + (110) + (20)}] 
K = 625) ~{(100) + (110) + (20) 1 


575 — 230 100 
= — + 
625 — 230 
K = 100 
= —— x 
395 
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Table 3. Sample points applied to classified image for the year 2005. 


Class Vegetation Build-up Water Total-User 
Vegetation 8 2 10 
Build-up Area 10 10 
Water 1 4 5 
Total-Producer 9 12 4 25 
_ (TS * TCS) — )i(col. total x Row total) 
Kappa coefficients(K ) = —————-—_____—Y—————_ * 100 


(TS)? — ¥ (col. total x Row total) 


g — [25+ 22) - {10 +9) + (10 +12) +S * 4H oo 
~~ (25)2—{(10* 9) + (10*12)+(5*4} ` 


— [(550) — {(90) + (120) + (20)}] 
K = 625) — (00 + (120) + 20} 0 


550-230 | 
= —C 
625 — 230 
=~ «100 
= * 
395 


K = 0.81 x 100 = 81% 


Table 4. Sample points applied to classified image for the year 1990. 


Class Vegetation Build-up Water Total-User 
Vegetation 9 1 10 
Build-up 9 1 10 
Water 1 4 5 
Total-Producer 10 10 5 25 
E (TS x FCS) — f (col. total « Row total) 
Kappa coefficients(K) = ——————-—___————_ * 100 


(TS)? — (col. total x Row total) 


ae — [25 * 22) ~ {(10 = 10) + (10 + 10) + (5531 oo 
~~ (25)2—{(10* 10) + (10 = 10) +(5*5)} 
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_ [(550) — {(100) + (100) + (25)}] 
K = (635) — {(100) + (100) + (25) 0 


550 — 225 ‘06 
= ———— * 
625 — 225 
K = 100 
= * 
400 


K = 0.8125 x 100 = 81.25% 


It determines that the classified image obtained using satellite image for the year 1990 
was 81.25%. The classified image for the year 2005 was 81% correct and the satellite image 
for the year 2019 was 87.34% correct. 


Discussions 


The results proved that most of urban sprawl occurred in Lahore and Gujranwala 
during 1990-2019. It was observed that mostly migrants rushed toward urban areas in search 
of education and employment. Lahore is known as an origin of education. Therefore, Lahore 


was a hotspot for knowledge seekers. 


People moved toward Lahore to get knowledge and stayed on rent with permanent 
residents. Soon, this trend become a business and the residents built new settlement to offer 
their houses on rent to migrants. In this way, about 85% houses were built without following 
any architectural by laws that devasted the complete face of the city. Most of migrants did not 
return to their native cities because they engaged themselves in some professional jobs and 
become permanent residents which cause to increase the urban area. 


Lahore and Gujranwala have limited health facilities in comparison to existing 


population. Mostly hospitals were built by British and a few were built after independence. 
Therefore, the increasing population is given cheap medical service due to lack of medical 
facilities. 


It is very difficult to apply law and order to ensure the security of masses with limited 
resources. Therefore, the crime rate had been increasing in Lahore and Gujranwala with 
passage of time since 1990. Transport is a remarkable problem of a mega city. It is very difficult 
to manage transportation of a big population. For this reason, many underpasses, overhead 
bridges, tunnels and metro lines were constructed which are serving within cities but traffic 
remains out of control. 
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On the other hand, Gujranwala is an industrial city. There were 550 industries in 
Gujranwala in 1995 that raised up to 2131 in 2019. These industries are famous for 


manufacturing goods of export quality which remained an attraction for labor and other 
professions to earn their livelihood which cause increase in urban size. 


A city faces a big burden of accommodation of migrants, hence a planned city built by 
senior architects is highly needed but unfortunately migrants built their houses by their own 
in Lahore and Gujranwala. Therefore, insecure and ill planned houses were built. 


Urban climates are getting polluted day by day by industrial emission and vehicle 
discharge. The statistics of Gujranwala describe that the amount of CO was 105 PPBV in 1995 
that increased up to 130 PPBV in 2018, CO2 increased from 363 PPM to 421 PPM. The 
number of vehicles were enhanced from 25730 in the year 2000 to 796751 in 2016. There were 
550 industries 1995 which increased up to 2131 in 2016. All these activities which are against 
the human health are a question mark to a step forward for a sustainable future. 


Conclusion. 

Urban sprawl is one of major problems globally for developing countries like Pakistan. 
It is addressed at every forum and require careful attention for immediate solutions. It is hard 
to draw city boundaries by physical visits therefore; remotely sensed datasets are highly reliable 
for identification of existing land use. But a limitation that we faced to fulfill the project 
successfully, was the coarse resolution of Landsat imagery. We executed field visits to 
determine the accuracy of results. One can get efficient results with improved satellite imagery 
of high resolution. 
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